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Abstract 



We present recent results on scalar light mesons based on Dalitz plot anal- 
yses of charm decays from Fermilab experiment E791. Low mass scalar mesons 
are found to have large contributions to the decays studied, D + — > K~ir + ir + 
and D + ,D+ — > tt^tt + it + . These results demonstrate the importance of charm 
decays as a new environment for the study of light meson physics. 



1 Introduction 

Here we present an overview of the results we obtained analysing the decays — > 
7r~7r + 7r + [I], D + —> 7T~7r + 7r + and D + — > K~7t + tt + ||, using data collected in 

1991/92 by Fermilab experiment E791 from 500GeV / ctt nucleon interactions. For 

details see [ref791]. 

To obtain these results, we had to introduce a new approach in Dalitz plot analysis 
in order to extract the mass and width of the scalar resonances by allowing them as 
floating parameters in the fit. We begin this paper presenting the general method, 
applied in the D + — > K~tt + tt + Dalitz-plot analysis, then we discuss the Df — > 
7r~7r + 7r + and D + — > studies using the same procedure. 



2 The D + i^-7r + 7T + Dalitz-plot Analysis 

From the original 2 x 10 1 events collected by E791, and after reconstruction and 
selection criteria, we obtained the k~7T + 7r + sample shown in Figure |l](a). The cross- 
hatched region contains the events selected for the Dalitz-plot analysis. There are 
15090 events in this sample, of which 6% are background. 

Figure |l](b) shows the Dalitz-plot for these events. The plot presents a rich struc- 
ture, where we can observe the clear bands from K (890)7r + , and an accumulation of 
events at the upper edge of the diagonal, due to heavier resonances. To study the res- 
onant substructure, we perform an unbinned maximun-likelihood fit to the data, with 



1 









N 3.5 












' (GcV/c 






j'i'il' 








f 




!(;!,'; 


1 1 lli'f'V'ii 1 


6 no 












4 DO 


J 


X 


0.5 







3 3.5 



■V (GeV/c7 



Figure 1: (a) The K 7r + 7r + invariant mass spectrum. The filled area is background; 
(b) Dalitz plot corresponding to the events in the dashed area of (a). 



probability distribution functions (PDF's) for both signal and background sources. 
In particular, for each candidate event, the signal PDF is written as the square of the 
total physical amplitude A and it is weighted for the acceptance across the Dalitz 
plot (obtained by Monte Carlo (MC)) and by the level of signal to background for 
each event, as given by the line shape of Figure |l|(a). The background PDF's (levels 
and shapes) are fixed for the Dalitz-plot fit, according to MC and data studies. 

We begin describing our first approach to fit the data, which represents the conven- 
tional Dalitz-plot analysis including the known Kir resonant amplitudes (A n , n>l), 
plus a constant non-resonant contribution. The resonance amplitudes used to describe 
the signal, Breit-Wigner parametrizations with Blatt-Weisskopf damping factors are 
described in reference [[|. 

Using this model with well-known resonances, we find contributions from the 
following channels: the non- resonant, responsible for more than 90% of the decay 
rate, followed by K*(1430)vr+, K*(892)vr+, K*(1680)vr+ and iT*(1430)7r + . The decay 
fractions and relative phases are available in reference 

To evaluate the fit quality, we compute a \ 2 from binned, two-dimensional distri- 
butions of data and decay model events. The \ 2 comes from the differences in the 
binned numbers of events between the data and the model (from a fast MC simu- 
lation). We obtain y^jv = 2.7 [y being the number of degrees of freedom), with a 
corresponding confidence level (CL) of 10" 11 . We thus conclude that a model with 
the known Kit resonances, plus a non-resonant amplitude, is not able to describe the 
D + — > K~7t + 7t + Dalitz plot satisfactorily. Thus, we are led to try an extra scalar 
resonance in our fit model. This second fit model, is constructed by the inclusion of 
an extra scalar state with unconstrained mass and width. For consistency, the mass 
and width of the other scalar state, the Kq(1430), are also free parameters of the fit. 
We adopt a better description for these scalar states by introducing gaussian-type 
form-factors || to take into account the finite size of the decaying mesons. 
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Figure 2: The n tt + ti + invariant mass spectrum. The dashed line represent the total 
background. Events used for the Dalitz analyses are in the hatched areas. 

Using this model, we obtain the values of 797 ± 19 ± 42 MeV/c 2 for the mass and 
410 ± 43 ± 85 MeV/c 2 for the width of the new scalar state (first error statistical, 
second error systematic), referred to here as the n. The values of mass and width 
obtained for the 7^(1430) are respectively 1459 ± 7 ± 6 MeV/c 2 and 175 ± 12 ± 12 
MeV/c 2 , appearing heavier and narrower than presented by the PDG 0. The decay 
fractions and relative phases with systematic errors, are given in reference ||. The 
K7r + state is now the dominant channel with decay fraction about 50%. 

3 The D+ 7T"7r + 7r + Results 

In Figure |]we show the 7r~7r + 7r + invariant mass distribution for the sample collected 
by E791 after reconstruction and selection criteria [II, 0. Besides combinatorial back- 
ground, reflections from the decays D + — > K~tt + tt + , D° —* K~n + (plus one extra 
7r + track) and — > rj'7r + , rj' — > p°(770)7 are all taken into account. The hatched 
regions in Figure show the samples used for the Dalitz-plot analyses. There are 1686 
and 937 candidate events for D + and respectively, with a signal to background 
ratio of about 2:1. The Dalitz plots for these events are shown in Figure ^, the axes 
corresponding to the two 7r _ 7r + invariant-masses squared. 

For the — > 7t~7t + tt + events in Figure 0(a), the signal amplitude includes all 
channels with well known dipion resonances |J: p°(770)7r + , / (980)7r + , / 2 (1270)7r + , 
/o(1370)7r + , p°(1450)7r + and the non-resonant, assumed constant across the Dalitz 
plot. 

The measured /o(980) standard Breit-Wigner parameters are m = 975 ± 3 
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Figure 3: (a) The —> n 7t + tt + Dalitz plot and (b) the D + — > n 7r + 7r + Dalitz 
plot. Since there are two identical pions, the plots are symmetrized. 

MeV/c 2 and T = 44±2±2 MeV/c 2 . The confidence level of the fit for D+ -> n-7i + n + 
is 35%. The decay fractions and relative phases, with systematic errors, are given in 
reference M. The /o(980)7r + state is the dominant channel with decay fraction about 
50%. 

In Figure |] we show the ii~7i + mass-squared projections for data (points) and 
model (solid lines, from fast-MC). 




Figure 4: (a) Su and S13 (m^) projections for Df — > 7r~7r + 7r + data (dots) and 
our best fit (solid), (b) S12 and S13 projections for D + — > 7r~7r + 7r + data (dots) and 
our best fit (solid) for models with a(500)7r + amplitude. In both figures the dashed 
distribution corresponds to the expected background level. 



4 The D + 7r"7r + 7r + Results 

In a first approach, we try to fit the D + — > 7r _ 7r + 7r + Dalitz plot of Figure 0(b) with 
the same amplitudes used for the Df — > tt~7t + tt + analysis. Using this model, the 
non-resonant, the p°(1450)7r + , and the p°(770)7r + amplitudes are found to dominate 
01 . However, this model does not describe the data satisfactorily, especially at low 
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7r~7r + mass squared 0. The \ 2 / v obtained from the binned Dalitz plot for this model 
is 1.6, with a CL less than 10 -5 . 

To investigate the possibility that another 7c~tt + resonance contributes to the 
D + — > 7r _ 7T + 7r + decay, we add an extra scalar resonance amplitude to the signal 
PDF, with mass and width as floating parameters in the fit. 

We find that this model improves our fit substantially. The mass and the width of 
the extra scalar state are found to be 478^23 ± 17 MeV/c 2 and 324±^ ± 21 MeV/c 2 , 
respectively. Refering to this state as the a (500), we observe that the cx(500)7r + 
channel produces the largest decay fraction 0; the non-resonant amplitude, which is 
dominant in the model without <j(500)7r + , drops substantially. The model with the 
sigma describes the data very well, as can be seen by the nn mass squared projection 
in Fig. |](b). The x 2 l v is now 0.9, with a corresponding confidence level of 91%. 



5 Conclusion 

From the data of the Fermilab E791 experiment, we studied the Dalitz plots of the 
decays D + — > K~it + tt + , — > 7r~7r + 7r + , and D + — ► 7r~7r + 7r + . In these three final 
states, the scalar intermediate resonances were found to give the largest contributions 
to the decay rates. We obtained strong evidence for the existence of the a (500) 
and k scalar mesons, measuring their masses and widths. We also obtained new 
measurements for masses and widths of the other scalars studied, /o(980), /o(1430) 
and if *(1430). 

The results presented here show the potential of D meson decays for the study of 
light meson spectroscopy, in particular in the scalar sector [[?]]. 
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